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We discuss the possibility to use a high energy neutrino beam from a muon storage ring to provide 
one way communication with a submerged submarine. Neutrino interactions produce muons which 
can be detected either, directly when they pass through the submarine or by their emission of 
Cerenkov light in sea water, which, in turn, can be exploited with sensitive photo detectors. Due to 
the very high neutrino flux from a muon storage ring, it is sufficient to mount either detection system 
directly onto the hull of the submersible. The achievable data transfer rates compare favorable with 
existing technologies and do allow for a communication at the usual speed and depth of submarines. 



The use of neutrino beams for communication is an old 
idea and has been put forth by several authors for vari- 
ous purposes, like for e.g interstellar or even intergalactic 
communication [l]-E]- Also the the use of neutrinos for 
communication with a submarine deep under the ocean 
has been previously considered [5 , 6 , however the conclu- 
sion was that this does not provide a feasible approach^. 
In this letter, we carefully reexamine the problem and 
will find that recent technological advances require to re- 
consider the earlier negative conclusion. 

Nuclear powered submarines offer, practically, unlim- 
ited submerged endurance; they are only tied to the sur- 
face by their need to communicate. Therefore, communi- 
cation at operational speed and depth is highly desirable. 
Currently, only radio transmission at extremely low fre- 
quency (elf) of < 100 Hz is able provide communication 
at speed and depth, elf data rates are very low, of 
order one bit per minute because of the very low band- 
width, the high noise levels and the difficulty to generate 
high-powered signals [8]. Instead radio transmission at 
frequencies of a few kHz (very low frequency, VLf) are 
used, providing data rates around 50bit/s [9 ; however, 
the sea water penetration of VLF is limited. This requires 
a wire antenna, floated close to the surface, which entails 
significant operational limitations. 

The basic concept for submarine neutrino communica- 
tion (SNC), derives from the fact that neutrinos can tra- 
verse an entire planet. Neutrinos can be sent from one 
point on the surface of the earth to a submarine, irrespec- 
tive of its location and depth. The feebleness of neutrino 
interactions implies a large detector and a very bright 
neutrino source. One of the currently most intense neu- 
trino beams is used in the MINOS experiment [10 , where 
a beam of neutrinos is sent from the Fermi National Ac- 
celerator Laboratory in Chicago to a mine in northern 
Minnesota over a distance of more than 700 km. Dur- 
ing the 2 years duration of the experiment, the neutrino 
beam produced 730 muons in 5 000 tons of detector. This 
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^ The use of hypothetical particles for this purpose is discussed in 
reference i7|. 



is equivalent to one neutrino event every 12 hours. Obvi- 
ously, an improvement of at least 6 orders of magnitude 
is required. 

Muon storage rings have been proposed as source of 
highly collimated neutrino beams [11 in order to allow 
precision measurements of the neutrino mixing param- 
eters [12]. In these facilities muons will be produced 
from pion decay and the pions are produced by proton 
irradiation of a target. Current designs for such a fa- 
cility assume lO^'^s-^ useful muon and 10^^ s ^ useful 
anti-muon decays with muon energies in the range from 
25—50 GeV. Such a facility also would constitute the first 
step towards a multi-TeV muon collider. The short life- 
time of the muon and the high proton beam power are 
the main technical challenges. Currently, several R&D 
experiments are underway, for a recent review see [13], 
in order to prove the feasibility of this concept. 

The muon neutrino flux from a storage ring with un- 
polarized muons of energy Eq is given by 

(1) 

where is the number of muon decays, P = 

yl — m^yEg , denotes the muon mass and L is the 
Cartesian distance from the storage ring. The beam di- 
vergence, 00 = tti^/Eq. These neutrinos will interact 
in sea water and produce muons in this process, with 
a cross section of da/dEj^ = ctq [Q -\- Q E'^/E'^^ , with 
cro = 1.583 • 10-^^2 m^, Q = 0.41 and Q = 0.08. The cross 
section for anti-neutrinos is obtained by interchanging Q 
and Q. The range of a muon, i?, with energy E^ in water 
can be parametrized [T4^ as R = 14.5 m+3.4 GeV~^ m 
in the range = 20 — 150 GeV. The resulting muon flux 
(j)^ per unit area is obtained from 

where p = 1.02 gcm~^ is the density of sea water, Na = 
6.3 • 10^^ is Avogadro's number and E_ is the smallest 
acceptable muon energy. The anti-muon flux is obtained 
by replacing the cross section and flux with the corre- 
sponding quantities for anti-neutrinos. The requirement 



that the muon direction is within less than a maximum 
angle Oj^ of the original neutrino direction, translates, by 
simple kinematics, into a lower bound Eq on the muon 
energy E^^ Eq = Ey/{1 + 2Ey9l). Also, we may set a 
minimum muon energy ^min every muon should have, in 
which case E- — maxj^^^*, ^min}- From equation [5] we 
obtain including both neutrinos and anti-neutrinos, with 



lOGeV, ^0 = 150GeV, AT^ 



10^ 



L = 10 000 km and = 10^ m^, which is the area in 
which a muon can be detected, a muon rate of = 
2.2 s~^; note, that (j)^ ex Eq. A neutrino source delivering 
muons at a rate of 10^^ s~^ with an energy of 150 GeV 
would require about 4MW in proton beam power and 
2.4 MW acceleration power, which for a 10% electrical 
efficiency translates into a total power consumption of 
roughly 65 MW. In order to aim the beam, it is nec- 
essary to point the long straight section of the muon 
storage ring towards the submarine. The storage ring is 
relatively large, but also quite lightweight and thus one 
solution could be to suspend the storage ring in water, 
either in a lake or close to shore, and use buoyant forces 
to aim it. Also, stronger magnets combined with shorter 
straight sections, can reduce the size of the ring substan- 
tially. The aiming of the beam represents a considerable 
engineering challenge but hardly can be considered an 
insurmountable obstacle. 

How many bits of information can be transmitted by 
one neutrino? SNC, in this respect, is very similar to 
deep space optical communications, where very few pho- 
tons need to carry the largest possible amount of infor- 
mation [15 . This is achieved by using pulse position 
modulation (ppm), where one unit of time is divided 
into M slots and we can freely chose in which of these 
slots we transmit the pulse. Of course, if enough pho- 
tons/neutrinos are available we can decide to transmit a 
number, P, of pulses per unit time. In the context of a 
muon storage ring, the number of slots is determined by 
the inverse of the duty factor. Duty factors of the order of 
10~^ seem feasible, therefore we will take M = 2^^. The 
optimum number of pulses will depend on the available 
muon event rate. In the context of information theory the 
above system corresponds, in the absence of backgrounds, 
to an M-ary Poisson erasure channel. The capacity of an 
information channel is the theoretically maximal rate at 
which information can be transmitted with an arbitrar- 
ily small error probability. In practice, a coding scheme 
is required to achieve good performance. However, re- 
search in the last decade or so has produced a number 
of practical coding schemes which actually can perform 
very close to capacity, see e.g. p^. The capacity, in our 
case, given as a function of P [17] is 
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where e is the erasure probability, which for a Poisson 
process is 1 — (1 — e^^^)^ with S being the number of 
events per unit time. For a given value of M and S we 
can determine the optimum number of pulses which max- 
imizes the capacity. For M = 2^^ the maximum capacity, 
Cm, can be parametrized as Cm = 6.61 • ^^"^^ bit /s. Using 
the previous result on (j)^ = 2.2 s~^ we obtain a capacity 
of ~ lObit/s even at the antipodes of the sender. 

Muons from cosmic ray interactions in the atmosphere 
and muons produced by neutrinos from the same cos- 
mic ray interactions are a source of irreducible back- 
ground. The muon flux at a water depth of 25 m is 

20s"^m"^sr"^ [14 , which for = 10^ m^ yields 
a background of 2 • 10^s~^sr~^. This flux, however, 
is entirely down going, i.e. there are no atmospheric 
muons from below the horizon. Thus, this background is 
problematic only as far as muons coming from the up- 
per hemisphere are erroneously reconstructed as com- 
ing from below. Neutrino telescopes, like antares or 
IceCube have demonstrated that the fraction of muons 
mis-reconstructed in this way is of the order 10~^ [18]. 
Moreover, for neutrinos stemming from a beam, we re- 
quire the resulting muons to lie within a cone with an 
opening angle of 20iy around the beam direction which 
yields another factor of 10 reduction. Thus, the resulting 
effective background rate, even close to the water surface 
is at most 2 • 10~^ s~^. The event rate from atmospheric 
neutrinos with energies above 10 GeV is a few events per 
day. Therefore, these backgrounds are at least 1000 times 
smaller than the signal and thus, can be safely neglected. 

The most straightforward approach to muon detection 
is to convert the submarine itself into a muon detector: a 
modern submarine is approximately a cylinder of 10 m di- 
ameter and 100 m length, which allows an effective muon 
detection area of 100 X 10 m^ = 10^ m^ We would 

use thin muon detector modules which can be used very 
much like wallpaper to cover the majority of the vessel's 
hull^. The muons would enter on one side of the subma- 
rine and leave it on the other side. The entry and exit 
points are measured and thus the the muon direction 
can be reconstructed quite precisely. Thin, flat and pos- 
sibly flexible, large area muon detectors exist in various 
forms like plastic scintillators, gas electron multipliers, 
see e.g. [19j, resistive plate chambers, etc. The require- 
ments in terms of timing and spatial resolution would be 



^ Note, that submarines have actually two hulls, a so called light 
hull, which defines the exterior, hydrodynamic shape and a pres- 
sure hull, which protects the crew and equipment from the hy- 
drostatic pressure. All detectors we consider, would be either 
mounted onto the light hull or in the space between the light 
and pressure hull. 
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well within the bounds of existing technologies and large 
area muon detectors are part of many high energy physics 
experiments. By exploiting local coincidence and coinci- 
dence across the vessel non-muon backgrounds like local 
radioactivity, either in the water or on-board the subma- 
rine, can be effectively rejected. Random coincidences 
from two or more muons striking opposite walls of the 
submarine can be effectively controlled by requiring the 
timing of the events to match the timing of an actual, 
single beam muon^. This approach of muon detection 
which we will refer to as direct^ has the advantage that 
the achievable and hence data transmission rate does 
not depend on environmental factors. The global distri- 
bution of data transmission rates for the direct detection 
scheme is shown in figure taking = 10^ m^. 

Muon detection in a transparent medium, like sea wa- 
ter, can be efficiently achieved by exploiting Cerenkov 
light, which is emitted along a cone with a constant 
opening angle, Oc] in water Oc = 42°. The number of 
photons emitted per unit length of the muon track in 
the wavelength range 400 — 550 nm, corresponding to 
the region of maximum transmission in sea water, is 
dN/dl 40 000m~^. The number of photons at the 
detector, which is a distance, r, away from the track 
is = eqexp[— r/ {sin OcX)] a/ {27rr) dN/dl, where A is 
the absorption length, a is the detector area and Cq is 
the quantum efficiency of the detector. Setting a trig- 
ger threshold of N = Nt, we can solve for the distance 
Vt up to which a muon will yield at least Nt photons. 
The effective muon detection area is then = r^ir. 
In order to allow track reconstruction, especially in the 
presence of backgrounds, it will be necessary to deploy 
a number, Q, of photon detectors, each with area a. In 
this case, the total number of signal photons, Sp will be 
Sp = NfQ and the total effective photon detection area, 
Ap will be aQ. 

Before we delve in detail into photonic backgrounds, it 
is useful to summarize the properties of the signal: 

1. Cerenkov light is emitted on a cone surrounding the 
muon track and in absence of scattering will travel 
in a straight line. Therefore, the relation between 
the arrival time of a photon at the photo-detector 
and its point of origin is unique. 

2. The arrival directions at each detector must sat- 
isfy the Cerenkov condition, i.e. the angle between 
the direction of muon and the photon must be the 
Cerenkov angle 0^ Furthermore, the arrival direc- 
tion, when traced back must intersect the muon 
track at one point. 



^ The rate of muon bundles, i.e. several muons originating from 
the same cosmic ray primary, is small enough not to affect our 
background estimates, see e.g. ,20^. 




FIG. 1: Global distribution of the achievable information ca- 
pacity. Panel (a) shows the Cerenkov detection scheme dur- 
ing daylight (noon), panel (b) shows the results for the di- 
rect detection scheme and panel (c) shows the results for the 
Cerenkov detection scheme during night. The black x de- 
notes the position of the sender. Land masses are dark gray 
and the turquoise area indicates missing data on the optical 
properties of sea water. The areas in olive indicate shallow 
water with a depth of 25 m or less. The maps are drawn in 
Hammer- Aitoff projection with a central meridian of 90° E. 

3. The mean number of photons at each detector is 
determined by the distance between the detector 
and the track. 

We will give a combinatorial argument for the approx- 
imate background rejection which can be obtained by 
exploiting these signatures. Assuming a timing resolu- 
tion of each detector of St we want each detector i to 
be no larger than dx = cSt. Conversely, each detector 
sees the track for no longer than 5t and thus we have a 
background suppression of St~^. Obviously, the detec- 
tor area, a, can be no larger than 5t^ tt. To keep the 
total number of photons detected constant it will be nec- 
essary to increase the number of detectors, Q, such that 
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Ap = Qa = Qc^dt'^TT stays constant. Also, individual 
detectors should be spaced no closer than Sx. At the 
same time, Sx sets the scale down to which the compo- 
nents of a track vector be determined. If we assume a 
linear array of photo-detectors with a length /p, then the 
angular resolution SO of the array is approximately given 
by tan SO = Sx/lp. Since we know that the muon must lie 
within SOjj of the beam direction, there can be, ne^ dif- 
ferent angles for a valid signal muon track uq = SO^jSO^ . 
The photon detectors survey an area with a maximal 
resolution of Sx^ thus there can be Ux different origins 
for the muon track Ux = A^^jSx^ . And finally, there can 
be different event starting times, nt = S^^. Altogether, 
there can he m = ng rix nt different valid muon tracks. 
For = 1 ns and Ip = 100 m, we obtain m ~ 10^^. 

Along each, of these m distinct, possible tracks there 
can be random fluctuation of stochastic backgrounds 
from either actual photons, radioactive decays or elec- 
tronic noise. Irrespective of their origin they all share 
the property of being random and thus are neither corre- 
lated in space nor time. The total background is given by 
B = AirbfAp, where bf is the background rate in units of 
s~^ m~^ sr~^. The track is seen at each individual photo- 
detector for only St, thus we need to be concerned only 
about background in the window St, i.e. we obtain an 
immediate background suppression of St~^ and the effec- 
tive background, = B St. We next require that the 
probability of the effective background B^ to fluctuate up 
to the signal Sp is smaller than {fm)~^, thus ensuring 
that on average there are nor more than 1// random co- 
incidence events per second happen in the entire detector 
and we use / = 1000. For large S and B, the probabil- 
ity can be approximated by a GauBian and we can solve 
for Be as a function of S and obtain B^ = S'^ /{2q^) with 
q = erfc~^ ^)] and where erfc~^ denotes the inverse 
of the complementary error function. 

So far, we have not made use of the photon arrival 
direction as discussed in point [2] Reconstruction of the 
photon arrival directions with a resolution of SOa will 
result in a background suppression, e^, of = SO'^. A 
simple optical system like a lens or mirror can translate 
the arrival direction of a photon into the position in the 
focal plane. Finally, the photon background can be as 
large as 

6<1.2ao»Sj(g)'(=5)s-.„-s.-.. (4) 

without causing a fake muon rate in excess of 1/f = 
10-^ s-^ For a given 6, we can set the trigger threshold 
Nt such that Sp is sufficiently large for equation [4] to hold. 
Nt, in turn, determines, together with A, how far out or 
at which maximal distance, r^, a muon can be seen above 
background and thus sets the size of A^. 

Next, we need to compare this rate to the various 
background light sources. The light produced by nat- 
ural radioactivity, mainly from ^^K, and the light from 



bio-luminescence is, at the depths we consider, negligible 
compared to the photons from the surface [21 . Nonethe- 
less, we include a random noise rate due to ^^K and 
bio-luminescence of 20kIIzcm~^, this corresponds to 20 
times the maximum noise reported in [21 . At a water 
depth, d, the photon flux, (j)p will be (j)p = (j)^ exp—d/X 
where A is the attenuation length and the photon 
flux at the surface. The photon flux from the sun in 
the zenith, in the wavelength range 400 — 550 nm, is 
3 ■ 10^^ s~^ m~^ sr~^. At night the surface pho- 
ton flux will be at least 25 000 times smaller. The larger 
the optical overburden the smaller Sp can be chosen and 
the larger the resulting A^ will be and vice versa. 

The optical properties of sea water play a crucial role 
for both the available target mass, and thus signal level 
as well as for the photon background. The attenua- 
tion length of sea water at a wavelength of 490 nm can 
be determined in situ by spectroscopy of the upwelling 
sun light. There are currently two satellites. Aqua and 
Terra, which carry on board the Moderate Resolution 
Imaging Spectroradiometers (modis) and provide nearly 
global coverage. This method has been validated by 
surface measurements, but, essentially, is sensitive only 
to a depth down to one attenuation length. There- 
fore, we will make the assumption that the attenuation 
length close to the surface is representative of the one 
at greater depth. We use a data set which has been av- 
eraged over 7 years and which has a spatial resolution 
of 2.5arcmin 4.6 km [22]. In order to estimate the 
available overburden we use global bathymetry data [23] 
with a horizontal resolution of 1 arcmin 2 km. For 
our numerical results we assume a diving depth of 300 m 
and a photon detection area Ap = 10 m^ with Cq = 0.5, 
St = Ins and SOa = 1°. We predict the expected infor- 
mation capacity: the results are shown for day time in 
figure [l^ and for night time in figure [Tj^. 

How such a system would be used operationally is dif- 
ficult to assess, since our knowledge of submarine opera- 
tions is insufficient to provide a in depth discussion. Ob- 
viously, the position of a submarine is secret and ideally 
not known to anyone outside the vessel. However, in or- 
der to send messages using a neutrino beam we need know 
the position of the submarine within the width of the 
beam of ~ 1 — 5 km. One solution, is to point the beam 
at prearranged times to prearranged locations within the 
patrol area of the submarine and the submarine will be 
at one of these points at the right time at least a cer- 
tain number of times per day. Another possibility, is to 
artificially increase the beam spot size by sweeping the 
beam across the suspected course of the submarine. The 
resulting data rate is low, but all that needs to transmit- 
ted is the request to sail to a prearranged location for 
the reception of high speed communication. Assuming a 
reduced data rate corresponding to the elf data rate of 1 
bit per minute, it follows from the expression for Cm that 
a muon rate of 3 • 10~^ s~^ is sufficient. This rate is more 
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than 7000 times smaller than the rate for a fully focused 
beam. The fully focused beam has an opening angle of 
6>o, corresponding to a beam spot area of about 40 km^ 
at distance of 10,000 km. Since we need only 1/7000 of 
the focused beam flux, we can sweep the beam over an 
area 7000 times larger than the beam spot area, which is 
7000 • 40 km^ = 280, 000 km^, corresponding to a sweep 
radius of o:^ 300 km. Note, that the sweep radius scales 
like ^0 and we have assumed direct detection only, opti- 
cal detection would yield significantly larger sweep radii. 
Finally, one could establish a neutrino mailbox^ consist- 
ing of a string moored to the seafloor, which has photo 
detectors attached, very similar to a single string of a 
conventional neutrino telescope like antares. The mail- 
box can be filled at very high data rate > 100 bit/s due to 
the deep sea location and resulting very low noise. The 
message is stored locally and can be retrieved by the sub- 
marine using short range communications. Each mailbox 
would be relatively cheap and thus a sizable number can 
be distributed throughout the patrol area. If one re- 
quires that the submarine is never further than 1 hour 
away from a mailbox and that the speed of a submarine 
is 20 knots, one obtains that each mailbox serves an area 
of 4, 000 km^. At 8 hour travel time an area similar to 
the one obtained from sweeping the beam can be covered 
by a single string and about 1400 strings would cover all 
oceans. Note, that the message size, due to the long in- 
tegration time of several hours combined with the high 
data rate, can reach many 100 kilobyte. 
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TABLE I: The fractional area of the world oceans in per 
cent for which the capacity given in the first column can be 
achieved. 

Table |T] summarizes our results succinctly. We have 
demonstrated that a neutrino beam from a muon storage 
ring can be detected by sensors mounted on the hull of 
a submarine. This in principle would allow to establish 
a one-way communication link at speed and depth with 
data rates in the range from 1 — 100 bit/s which improves 
current ELF data rates by 1-3 orders of magnitude and is 
similar to data rates offered by vlf. 
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